Studying the Large Scale Structure and interstellar Medium of Galaxies During the Epochs of Peak Cosmic Star formation and Reionization With Infrared Fine Structure Lines

Infrared (IR) fine-structure (FS) lines from trace metals in the interstellar medium (ISM) of galaxies are valuable diagnostics of the physical conditions in a broad range of astrophysical environments, such gas irradiated by stellar far-ultraviolet (FUV) photons or X-rays from accreting supermassive black holes, called active galactic nuclei (AGN). The transparency of these lines to dust and their high escape fractions into the intergalactic medium (IGM) render them as useful probes to study the epochs of peak cosmic star formation (SF) and Reionization. Chapter 1 of this thesis is a study of the ISM of the Cloverleaf quasar. Observations of IR FS lines from singly ionized carbon and neutral oxygen have allowed us to assess the physical conditions—parametrized by their gas density and the impingent FUV flux—prevalent in atomic gas heated by stellar FUV photons. We find that UV heating from local SF is not sufficient to explain the measured FS and molecular luminosities, and suggest that X-ray heating from the AGN is required to simultaneously explain both sets of data. The general picture of the Cloverleaf ISM that emerges from our composite model is one where the [CII] and [OI]63 line emission is produced primarily within PDRs and HII regions of a 1.3-kpc wide starburst, which is embedded in a denser XDR component that is the dominant source of heating for the CO gas. The fact that the star-forming PDR and HII region gas is co-spatial with the XDR—and within ∼ 650 pc of the accreting black hole—provides strong evidence that SF is ongoing while immersed in a strong X-ray radiation field provided by the nearby AGN. This finding has implications for the co-evolution of supermassive black holes and their host galaxies. The work in this chapter will be submitted for first-author publication imminently. In Chapter 2, we explore the possibility of studying the redshifted far-IR fine-structure line emission using the three-dimensional (3-D) power spectra obtained with an imaging spectrometer. The intensity mapping approach measures the spatio-spectral fluctuations due to line emission from all galaxies, including those below the individual detection threshold. The technique provides 3-D measurements of galaxy clustering and moments of the galaxy luminosity function. Furthermore, the linear portion of the power spectrum can be used to measure the total line emission intensity including all sources through cosmic time with redshift information naturally encoded. Total line emission, when compared to the total star formation activity and/or other line intensities reveals evolution of the interstellar conditions of galaxies in aggregate. As a case study, we consider measurement of [CII] autocorrelation in the 0.5 \u3c z \u3c 1.5 epoch, where interloper lines are minimized, using far-IR/submm balloon-borne and future space-borne instruments with moderate and high sensitivity, respectively. In this context, we compare the intensity mapping approach to blind galaxy surveys based on individual detections. We find that intensity mapping is nearly always the best way to obtain the total line emission because blind, wide-field galaxy surveys lack sufficient depth and deep pencil beams do not observe enough galaxies in the requisite luminosity and redshift bins. Also, intensity mapping is often the most efficient way to measure the power spectrum shape, depending on the details of the luminosity function and the telescope aperture. The work in this chapter has been published in Uzgil et al. (2014). In the final Chapter, we consider the extension of intensity mapping experiments targeting IR FS lines to the late stages of the Epoch of Reionization (EoR), at z ∼ 7. Intensity mapping experiments of emission lines from the ISM of galaxies are highly complementary to experiments that are aiming to detect the 21 cm power spectrum during the same epoch, as the former is a direct probe of the sources of Reionization, and the latter is a probe of the effect of those sources on the surrounding IGM. Since current and planned observations are limited by cosmic variance at the bright end of the galaxy luminosity function, and will not be able to detect the faintest galaxies responsible for a significant fraction of the ionizing photon supply during EoR, intensity mapping is an appealing approach to study the nature and evolution of galaxies during this stage in the history of the Universe. Again, the utility of FS lines as ISM diagnostics, combined with the ability of intensity mapping to measure redshift-evolution in mean intensity of individual lines or the evolution of line ratios (constructed from multiple cross-power spectra), presents a unique and tantalizing opportunity to directly observe changes in properties of interstellar medium (such as hardness of the ionizing spectrum in galaxies and metallicity) that are important to galaxy evolution studies

Constraining the ISM Properties of the Cloverleaf Quasar Host Galaxy with Herschel Spectroscopy

We present Herschel observations of the far-infrared (FIR) fine-structure (FS) lines [C II]158 μm, [O I]63 μm, [O III]52 μm, and [Si II]35 μm in the z = 2.56 Cloverleaf quasar, and combine them with published data in an analysis of the dense interstellar medium (ISM) in this system. Observed [C II]158 μm, [O I]63 μm, and FIR continuum flux ratios are reproduced with photodissociation region (PDR) models characterized by moderate far-ultraviolet (FUV) radiation fields with G_0 = 0.3–1 × 10^3 and atomic gas densities n_H = 3–5 × 10^3 cm^(−3), depending on contributions to [C II]158 μm from ionized gas. We assess the contribution to the [C II]158 μm flux from an active galactic nucleus (AGN) narrow line region (NLR) using ground-based measurements of the [N II]122 μm transition, finding that the NLR can contribute at most 20%–30% of the observed [C II]158 μm flux. The PDR density and far-UV radiation fields inferred from the atomic lines are not consistent with the CO emission, indicating that the molecular gas excitation is not solely provided via UV heating from local star formation (SF), but requires an additional heating source. X-ray heating from the AGN is explored, and we find that X-ray-dominated region (XDR) models, in combination with PDR models, can match the CO cooling without overproducing the observed FS line emission. While this XDR/PDR solution is favored given the evidence for both X-rays and SF in the Cloverleaf, we also investigate alternatives for the warm molecular gas, finding that either mechanical heating via low-velocity shocks or an enhanced cosmic-ray ionization rate may also contribute. Finally, we include upper limits on two other measurements attempted in the Herschel program: [C II]158 μm in FSC 10214 and [O I]63 μm in APM 08279+5255

Constraining the ISM Properties of the Cloverleaf Quasar Host Galaxy with Herschel Spectroscopy

We present Herschel observations of the far-infrared (FIR) fine-structure (FS) lines [C II]158 μm, [O I]63 μm, [O III]52 μm, and [Si II]35 μm in the z = 2.56 Cloverleaf quasar, and combine them with published data in an analysis of the dense interstellar medium (ISM) in this system. Observed [C II]158 μm, [O I]63 μm, and FIR continuum flux ratios are reproduced with photodissociation region (PDR) models characterized by moderate far-ultraviolet (FUV) radiation fields with G_0 = 0.3–1 × 10^3 and atomic gas densities n_H = 3–5 × 10^3 cm^(−3), depending on contributions to [C II]158 μm from ionized gas. We assess the contribution to the [C II]158 μm flux from an active galactic nucleus (AGN) narrow line region (NLR) using ground-based measurements of the [N II]122 μm transition, finding that the NLR can contribute at most 20%–30% of the observed [C II]158 μm flux. The PDR density and far-UV radiation fields inferred from the atomic lines are not consistent with the CO emission, indicating that the molecular gas excitation is not solely provided via UV heating from local star formation (SF), but requires an additional heating source. X-ray heating from the AGN is explored, and we find that X-ray-dominated region (XDR) models, in combination with PDR models, can match the CO cooling without overproducing the observed FS line emission. While this XDR/PDR solution is favored given the evidence for both X-rays and SF in the Cloverleaf, we also investigate alternatives for the warm molecular gas, finding that either mechanical heating via low-velocity shocks or an enhanced cosmic-ray ionization rate may also contribute. Finally, we include upper limits on two other measurements attempted in the Herschel program: [C II]158 μm in FSC 10214 and [O I]63 μm in APM 08279+5255

A Foreground Masking Strategy for [CII] Intensity Mapping Experiments Using Galaxies Selected by Stellar Mass and Redshift

Intensity mapping provides a unique means to probe the epoch of reionization (EoR), when the neutral intergalactic medium was ionized by the energetic photons emitted from the first galaxies. The [CII] 158$\mu$m fine-structure line is typically one of the brightest emission lines of star-forming galaxies and thus a promising tracer of the global EoR star-formation activity. However, [CII] intensity maps at $6 \lesssim z \lesssim 8$ are contaminated by interloping CO rotational line emission ($3 \leq J_{\rm upp} \leq 6$) from lower-redshift galaxies. Here we present a strategy to remove the foreground contamination in upcoming [CII] intensity mapping experiments, guided by a model of CO emission from foreground galaxies. The model is based on empirical measurements of the mean and scatter of the total infrared luminosities of galaxies at $z 10^{8}\,\rm M_{\rm \odot}$ selected in $K$-band from the COSMOS/UltraVISTA survey, which can be converted to CO line strengths. For a mock field of the Tomographic Ionized-carbon Mapping Experiment (TIME), we find that masking out the "voxels" (spectral-spatial elements) containing foreground galaxies identified using an optimized CO flux threshold results in a $z$-dependent criterion $m^{\rm AB}_{\rm K} \lesssim 22$ (or $M_{*} \gtrsim 10^{9} \,\rm M_{\rm \odot}$) at $z < 1$ and makes a [CII]/CO$_{\rm tot}$ power ratio of $\gtrsim 10$ at $k=0.1$ $h$/Mpc achievable, at the cost of a moderate $\lesssim 8\%$ loss of total survey volume.Comment: 14 figures, 4 tables, re-submitted to ApJ after addressing reviewer's comments. Comments welcom

No Evidence for Enhanced [O iii] 88 μm Emission in a z ∼ 6 Quasar Compared to Its Companion Starbursting Galaxy

We present Atacama Large Millimeter/submillimeter Array band 8 observations of the [O III] 88 mu m line and the underlying thermal infrared continuum emission in the z = 6.08 quasar CFHQS J2100-1715 and its dust-obscured starburst companion galaxy (projected distance: similar to 60 kpc). Each galaxy hosts dust-obscured star formation at rates >100 M-circle dot yr(-1), but only the quasar shows evidence for an accreting 10(9) M-circle dot black hole. Therefore we can compare the properties of the interstellar medium in distinct galactic environments in two physically associated objects, similar to 1 Gyr after the big bang. Bright [O III] 88 mu m emission from ionized gas is detected in both systems; the positions and linewidths are consistent with earlier [C II] measurements, indicating that both lines trace the same gravitational potential on galactic scales. The [O III] 88 mu m/far-infrared (FIR) luminosity ratios in both sources fall in the upper range observed in local luminous infrared galaxies of similar dust temperature, although the ratio of the quasar is smaller than in the companion. This suggests that gas ionization by the quasar (expected to lead to strong optical [0 III] 5008 angstrom emission) does not dominantly determine the quasar's FIR [O III] 88 mu m luminosity. Both the inferred number of photons needed for the creation of O++ and the typical line ratios can be accounted for without invoking extreme (top-heavy) stellar initial mass functions in the starbursts of both sources.ERC [740246]; National Science Foundation [AST-1614213]; National Science Foundation of China [11721303]This item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at [email protected]

Detector modules and spectrometers for the TIME-Pilot [CII] intensity mapping experiment

This proceeding presents the current TIME-Pilot instrument design and status with a focus on the close-packed modular detector arrays and spectrometers. Results of laboratory tests with prototype detectors and spectrometers are discussed. TIME-Pilot is a new mm-wavelength grating spectrometer array under development that will study the Epoch of Reionization (the period of time when the first stars and galaxies ionized the intergalactic medium) by mapping the fluctuations of the redshifted 157:7 μm emission line of singly ionized carbon ([CII]) from redshift z ~ 5:2 to 8:5. As a tracer of star formation, the [CII] power spectrum can provide information on the sources driving reionization and complements 21 cm data (which traces neutral hydrogen in the intergalactic medium). Intensity mapping provides a measure of the mean [CII] intensity without the need to resolve and detect faint sources individually. We plan to target a 1 degree by 0.35 arcminute field on the sky and a spectral range of 199-305 GHz, producing a spatial-spectral slab which is 140 Mpc by 0.9 Mpc on-end and 1230 Mpc in the redshift direction. With careful removal of intermediate-redshift CO sources, we anticipate a detection of the halo-halo clustering term in the [CII] power spectrum consistent with current models for star formation history in 240 hours on the JCMT. TIME-Pilot will use two stacks of 16 parallel-plate waveguide spectrometers (one stack per polarization) with a resolving power R ~ 100 and a spectral range of 183 to 326 GHz. The range is divided into 60 spectral channels, of which 16 at the band edges on each spectrometer serve as atmospheric monitors. The diffraction gratings are curved to produce a compact instrument, each focusing the diffracted light onto an output arc sampled by the 60 bolometers. The bolometers are built in buttable dies of 8 (low freqeuency) or 12 (high frequency) spectral channels by 8 spatial channels and are mated to the spectrometer stacks. Each detector consists of a gold micro-mesh absorber and a titanium transition edge sensor (TES). The detectors (1920 total) are designed to operate from a 250 mK base temperature in an existing cryostat with a photon-noise-dominated NEP of ~2 * 10^(-17) WHz^(-1-2). A set of flexible superconducting cables connect the detectors to a time-domain multiplexing SQUID readout system

The evolution of the cosmic molecular gas density

One of the last missing pieces in the puzzle of galaxy formation and evolution through cosmic history is a detailed picture of the role of the cold gas supply in the star-formation process. Cold gas is the fuel for star formation, and thus regulates the buildup of stellar mass, both through the amount of material present through a galaxy's gas mass fraction, and through the efficiency at which it is converted to stars. Over the last decade, important progress has been made in understanding the relative importance of these two factors along with the role of feedback, and the first measurements of the volume density of cold gas out to redshift 4, (the "cold gas history of the Universe") has been obtained. To match the precision of measurements of the star formation and black-hole accretion histories over the coming decades, a two orders of magnitude improvement in molecular line survey speeds is required compared to what is possible with current facilities. Possible pathways towards such large gains include significant upgrades to current facilities like ALMA by 2030 (and beyond), and eventually the construction of a new generation of radio-to-millimeter wavelength facilities, such as the next generation Very Large Array (ngVLA) concept.Comment: 7 pages, 2 figures, Science White paper submitted to Astro2020 Decadal Surve

Probing Cosmic Reionization and Molecular Gas Growth with TIME

Line intensity mapping (LIM) provides a unique and powerful means to probe cosmic structures by measuring the aggregate line emission from all galaxies across redshift. The method is complementary to conventional galaxy redshift surveys that are object-based and demand exquisite point-source sensitivity. The Tomographic Ionized-carbon Mapping Experiment (TIME) will measure the star formation rate (SFR) during cosmic reionization by observing the redshifted [CII] 158$\mu$m line ($6 \lesssim z \lesssim 9$) in the LIM regime. TIME will simultaneously study the abundance of molecular gas during the era of peak star formation by observing the rotational CO lines emitted by galaxies at $0.5 \lesssim z \lesssim 2$. We present the modeling framework that predicts the constraining power of TIME on a number of observables, including the line luminosity function, and the auto- and cross-correlation power spectra, including synergies with external galaxy tracers. Based on an optimized survey strategy and fiducial model parameters informed by existing observations, we forecast constraints on physical quantities relevant to reionization and galaxy evolution, such as the escape fraction of ionizing photons during reionization, the faint-end slope of the galaxy luminosity function at high redshift, and the cosmic molecular gas density at cosmic noon. We discuss how these constraints can advance our understanding of cosmological galaxy evolution at the two distinct cosmic epochs for TIME, starting in 2021, and how they could be improved in future phases of the experiment.Comment: 30 pages, 18 figures, accepted for publication in Ap

The ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Evolution of the Molecular Gas in CO-selected Galaxies

We analyze the interstellar medium properties of a sample of 16 bright CO line emitting galaxies identified in the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Large Program. This CO−selected galaxy sample is complemented by two additional CO line emitters in the UDF that are identified based on their MultiUnit Spectroscopic Explorer (MUSE) optical spectroscopic redshifts. The ASPECS CO−selected galaxies cover a larger range of star formation rates (SFRs) and stellar masses compared to literature CO emitting galaxies at z > 1 for which scaling relations have been established previously. Most of ASPECS CO-selected galaxies follow these established relations in terms of gas depletion timescales and gas fractions as a function of redshift, as well as the SFR–stellar mass relation (“galaxy main sequence”). However, we find that ∼30% of the galaxies (5 out of 16) are offset from the galaxy main sequence at their respective redshift, with ∼12% (2 out of 16) falling below this relationship. Some CO-rich galaxies exhibit low SFRs, and yet show substantial molecular gas reservoirs, yielding long gas depletion timescales. Capitalizing on the well-defined cosmic volume probed by our observations, we measure the contribution of galaxies above, below, and on the galaxy main sequence to the total cosmic molecular gas density at different lookback times. We conclude that main-sequence galaxies are the largest contributors to the molecular gas density at any redshift probed by our observations (z ∼ 1−3). The respective contribution by starburst galaxies above the main sequence decreases from z ∼ 2.5 to z ∼ 1, whereas we find tentative evidence for an increased contribution to the cosmic molecular gas density from the passive galaxies below the main sequenc